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(Dated: February 8, 2008) 

Elliptic flow holds much promise for studying the early-time thermalization attained in ultra- 
relativistic nuclear collisions. Flow measurements also provide a means of distinguishing between 
hydrodynamic models and calculations which approach the low density (dilute gas) limit. Among 
the effects that can complicate the interpretation of elliptic flow measurements are azimuthal cor- 
relations that are unrelated to the reaction plane (non-flow correlations). Using data for Au + Au 
collisions at ^/snn = 130 GeV from the STAR TPC, it is found that four-particle correlation anal- 
yses can reliably separate flow and non-flow correlation signals. The latter account for on average 
about 15% of the observed second-harmonic azimuthal correlation, with the largest relative contri- 
bution for the most peripheral and the most central collisions. The results are also corrected for the 
effect of flow variations within centrality bins. This effect is negligible for all but the most central 
bin, where the correction to the elliptic flow is about a factor of two. A simple new method for two- 
particle flow analysis based on scalar products is described. An analysis based on the distribution 
of the magnitude of the flow vector is also described. 



PACS numbers: 25.75.Ld 



I. INTRODUCTION 



In non-central heavy- ion collisions, the initial spatial 
deformation due to geometry and the pressure developed 
early in the collision causes azimuthal momentum-space 
anisotropy, which is correlated with the reaction plane 
[jl], p[ |[ [IJ. Measurements of this correlation, known 



as anisotropic transverse flow, provide insight into the 
evolution of the early stage of a relativistic heavy-ion 
collision || . Elliptic flow is characterized by the second 
harmonic coefficient vi of an azimuthal Fourier decom- 
position of the momentum distribution (6| [7| , and has 
been observed and extensively studied in nuclear colli- 
sions from sub-relativistic energies on up to RHIC. At 
top AGS and SPS energies, elliptic flow is inferred to be a 
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relative enhancement of emission in the plane of the reac- 
tion. Elliptic flow is developed mostly in the first several 
fm/c (of the order of the size of nuclei) after the colli- 
sion and thus provides information about the early-time 
thermalization achieved in the collisions ||. Generally 
speaking, large values of flow are considered signatures 
of hydrodynamic behavior r, |To|, although an alter- 
native approach [y| |i~4~15|, |l(f"is also argued to be 
consistent with the large elliptic flow for pions and pro- 
tons at RHIC H^j. Models in which the colliding nuclei 
resemble interacting volumes of dilute gas — the low den- 
sity limit jl8| (LDL) — represent the limit of mean free 
path that is the opposite of hydrodynamics. It remains 
unclear to what extent the LDL picture can describe the 
data at RHIC, and valuable insights can be gained from 
mapping out the conditions under which hydrodynamic 
and LDL calculations can reproduce the measured ellip- 
tic flow. 

Anisotropic flow refers to correlations in particle emis- 
sion with respect to the reaction plane. The reac- 
tion plane orientation is not known in experiment, and 
anisotropic flow is usually reconstructed from the two- 
particle azimuthal correlations. But there are several 
possible sources of azimuthal correlations that are un- 
related to the reaction plane — examples include corre- 
lations caused by resonance decays, (mini)jets, strings, 
quantum statistics effects, final state interactions (par- 
ticularly Coulomb effects), momentum conservation, etc. 
The present study does not distinguish between the var- 
ious effects in this overall category, but classifies their 
combined effect as "non-flow" correlations. 

Conventional flow analyses are equivalent to averag- 
ing over correlation observables constructed from pairs 
of particles. When such analyses are applied to relativis- 
tic nuclear collisions where particle multiplicities can be 
as high as a few thousand, the possible new informa- 
tion contained in multiplets higher than pairs remains 
untapped. A previous study of high-order flow effects 
focused on measuring the extent to which all fragments 
contribute to the observed flow signal [HJ , and amounted 
to an indirect means of separating flow and non-flow cor- 
relations. Given that flow analyses based on pair cor- 
relations are sensitive to both flow and non-flow effects, 
the present work investigates correlation observables con- 
structed from particle quadruplets. The cumulant for- 
malism removes the lower-order correlations which are 
present among any set of four particles, leaving only the 
effect from the so-called "pure" quadruplet correlation. 
The simplest cumulant approach, in terms of both con- 
cept and implementation, partitions observed events into 
four subevents. In the present study, the four-subevent 
approach is demonstrated, but our main focus is on a 
more elaborate cumulant method, developed by Borgh- 
ini, Dinh and Ollitrault [^0|, pi] . There are indications 
that non-flow effects contribute at a negligible level to 
the four-particle cumulant correlation p(i[ pi] , making it 
unnecessary to continue to even higher orders for the pur- 
pose of separating the flow and non-flow signals. This 



observation is confirmed by our Monte-Carlo simulations 
In this analysis the observed multiplicity of charged 
particles within the detector acceptance is used to char- 
acterize centrality. This leads to some fluctuations of the 
impact parameter and, correspondingly, of the elliptic 
flow within each centrality bin, especially in the bin of 
highest multiplicity. In the present study, a correction is 
applied to reduce a possible bias in the measurements of 
the mean elliptic flow due to impact parameter fluctua- 
tions in the centrality bins to an insignificant level. 

The present study begins with a review of the standard 
pair correlation method, and provides details concerning 
the approach adopted in earlier STAR publications [Tt]] 
for treating non-flow correlations. A new method of pair 
flow analysis using the scalar product of flow vectors also 
is introduced. In the conventional method, a flow coeffi- 
cient is calculated by the mean cosine of the difference in 
angle of two flow vectors. In the scalar product method, 
this quantity is weighted by the lengths of the vectors. 
The new method offers advantages, and is also simple to 
apply. Also, an analysis in terms of the distribution of 
the magnitude of the flow vector is discussed. 
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FIG. 1: The azimuthal angle distribution of tracks from 
minimum bias events. Dips are due to the reduced efficiency 
at sector boundaries of STAR TPC. 

Measurements presented in this paper are based on 
Au+Au data at ^/sWn = 130 GeV recorded by STAR 
(Solcnoidal Tracker At RHIC) during the summer of 
2000. A detailed description of the detector in its year- 
one configuration can be found elsewhere [ p2| . The main 
feature of the STAR TPC relevant to this analysis is 
its full azimuthal coverage (see Fig. |l|). The analysis 
is based on 170000 events corresponding to a minimum 
bias trigger. Events with a primary vertex beyond 1 cm 
radially from the beam or 75 cm longitudinally from the 
center of the Time Projection Chamber (TPC) were ex- 
cluded. Within the selected events, tracks were used for 
the estimation of the flow vector if all five of the follow- 
ing conditions were satisfied: they passed within 2 cm of 
the primary vertex, they had at least 15 space points in 
the TPC, the ratio of the number of space points to the 
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expected maximum number of space points was greater 
than 0.52, pseudorapidity \r)\ < 1.3, and transverse mo- 
mentum 0.1 < pt < 2.0 GeV/c. Particles over a wider 
range in r\ and pt were correlated with this flow vector as 
shown in the graphs below. Centrality is characterized in 
eight bins of charged particle multiplicity, n c h, divided by 
the maximum observed charged multiplicity, n max , with 
a more stringent cut \rf\ < 0.75 imposed only for this 
centrality determination. The above cuts are essentially 
the same as used in the previous STAR studies of elliptic 
flow @0. 



II. TWO-PARTICLE CORRELATION 
METHODS 

Anisotropic transverse flow manifests itself in the dis- 
tribution of <j)' = (f>—^t, where <f> is the measured azimuth 
for a track in detector coordinates, and W is the azimuth 
of the estimated reaction plane in that event. The ob- 
served anisotropics are described by a Fourier expansion, 



dN/d4>' oc 1 + 2vi_ b s cos </>' + 2v2,obs cos 2<f>' + 



(1) 



Each measurable harmonic can yield an independent es- 
timate of the event reaction plane via the event flow 
vector Q n : 



Qn COS n^n — COS TKpj , 

i 

Q n sin n\&„ = 2J sm n 4>i 7 



(2) 



where the sums extend over all particles in a given event. 
The observed values of v n b s corrected for the reaction 
plane resolution yield v n ^. Below we will also use the 
representation of the flow vector as a complex number 
with real and imaginary parts equal to x and y compo- 
nents defined in Eq. (Q): 



Qn = '' 



(3) 



where u n ,i = e m< " is a unit vector associated with the 
i-th particle; its complex conjugate is denoted by u ni . 



A. Correlation between flow angles from different 
subevents. Estimate of non-flow effects. 

In order to report anisotropic flow measurements in a 
detector-independent form, it is customary to divide each 
event into two subevents and determine the resolution 
of the event plane by correlating the Q n vector for the 
subevents f| gjj]. In order to estimate the contribution 
from different non-flow effects one can use different ways 
of partitioning the entire event into two subevents. The 
partition according to particle charge should be more af- 
fected by resonance decay effects because the decay prod- 
ucts of neutral resonances have opposite charge. The 



partition using two (pseudo) rapidity regions (better sep- 
arated by Ay > 0.1) should greatly suppress the con- 
tribution from quantum statistics effects and Coulomb 
(final state) interactions. 

Another important observation for the estimate of the 
non-flow effects is their dependence on centrality. The 
correlation between two subevent flow angles is 



cos(2(^ a) -^))) 



/ Z^=l ^ . 2^3 = 1 "j 

V Msub yMsub 

M su f)M su b 
= M h {UiU i } 

lvl sub 

cx M sub (v 2 2 +S 2 ), 



(4) 



where M su b is the multiplicity of a sub-event, and 82 
denotes the non-flow contribution to two-particle corre- 
lations. For correlations due to small clusters, which are 
believed responsible for the dominant non-flow correla- 
tions |2(J , the strength of the correlation should scale in 
inverse proportion to the total multiplicity. Since the 
subevent multiplicity is proportional to the total multi- 
plicity, we can define S 2 to be the multiplicity indepen- 
dent non-flow effect: 82 — 82/M su b- Collecting terms, we 
arrive at 



cos(2(* 2 a) - K MsuhV 2 + 62 _ 



(•5) 



What is important is that the non-flow contribution to 
(cos(2(5 , 2 !1 ' l ~^ , 2 b ' ) ))) 1S approximately independent of cen- 
trality. The typical shape of (cos(2(*^ q) -* 2 b) ))) for flow 
(see, for example, Fig. g) is peaked at mid-central events 
due to the fact that for peripheral collisions, Af su b is 
small, and for central events, V2 is small. In the previous 
estimates J^, [lTj of the systematic errors, we have set the 
quantity 82 = 0.05. The justification for this value was 
the observation of similar correlations for the first and 
higher harmonics (we have investigated up to the sixth 
harmonic). One could expect the non-flow contribution 
to be of similar order of magnitude for all these harmon- 
ics, and HIJING J^4| simulations support this conclusion. 
Given the value 82 = 0.05, one simply estimates the con- 
tribution from non-flow effects to the measurement of v 2 
from the plot of (cos(2(* 2 a) - l^))) using Eq. (§). 

Figured shows the event plane correlation between two 
subevents, for each of three different subevent partitions. 
In central events, it is seen that the correlation is stronger 
in the case of subevents with opposite sign of charge com- 
pared to subevents partitioned randomly. This pattern 
might be due to resonance decays to two particles with 
opposite charge. The spread of the results for different 
subevent partitions is about 0.05, which is in accord with 
the number used for the estimates of the systematic er- 
rors. 

The event plane resolution for full events is defined as 

(C0s(7l(*„ leas „ re -* trU e))>, hi which * measure and * true 

are azimuthal angles for the measured reaction plane and 
the "true" reaction plane, respectively. The resolution 



with p t weighting (see Section II B) can reach as high as 



5 



■Q CM 



ra cm 

> 

o 
o 

V 




0.09 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



n ch /n max 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



n ch /n max 



FIG. 2: Correlation between the event plane angles deter- 
mined from pairs of subevents partitioned randomly (circles), 
partitioned with opposite sign of pseudorapidity (squares) and 
partitioned with opposite sign of charge (crosses). The cor- 
relation is plotted as a function of centrality, namely charged 
particle multiplicity n c h divided by the maximum observed 
charged multiplicity, n max . 



FIG. 4: Elliptic flow signal V2 as a function of centrality, 
from study of the correlation between particle pairs consisting 
of randomly chosen particles (circles), particles with opposite 
sign of charge (crosses), particles with the same sign of charge 
(triangles) and particles with opposite sign of pseudorapidity 
(squares). 



0.8, as shown in Fig. |l The v-i as a function of cen- 
trality is shown in Fig.^j, using different prescriptions to 
partition the particles into subevents. Again, partition- 
ing into subevents with opposite sign of charge yields the 
highest elliptic flow signal, presumably because of neutral 
resonance (p°, etc.) decay. 
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FIG. 3: The event plane resolution for full events as a func- 
tion of centrality, using randomly partitioned subevents with 
(circles) and without (triangles) pt weight. 



B. Weighting 

If Eq. (|^) is generalized to the form Q n = J2i w i u ij 
where the Wi are weights adjusted to optimize the event 
plane resolution || ^5|, then it, should be replaced by 
WiUi for all equations in this paper, and M should be 
replaced by X^i^i throughout Sec. O, and Sec. IV B. 



The best weight Wi(r),p t ) is V2(r),Pt) itself Q. In prac- 
tice, since we know that vi is approximately proportional 
to p t up to about 2 GeV/c, it is convenient to use pt as 
the weight. It is found that pt weighting can reduce the 
statistical error significantly, as demonstrated in Fig. 



C. Scalar product flow analysis 

In a new scalar product method p6|] , each event is par- 
titioned into two subevents, labeled by the superscripts 
a and b. The correlation between two subevents is 



(Q°Q b n ) = (v 2 n M a M» 



(6) 



where M a and M b are the multiplicities for subevents 
a and b, respectively. The vectors and Q b n are con- 
structed for the appropriate subevent as per Eq. 

Given the above, the flow relative to the true reaction 
plane can be readily calculated from unit momentum vec- 
tors u nt i(r),pt) of the analyzed tracks by using Eq. (||) for 
the particle relative to the 2M other particles, and then 
dividing by the square root of Eq. (^|) for the subevents. 
This gives 



(QnU* nii (r],p t )) 



(J) 



Auto-correlations are removed by subtracting particle i in 
the calculation of Q n when taking the scalar product with 
u n> i. This method weights events with the magnitude of 
the Q n vector, and if Q n is replaced by its unit vector, 
the above reduces to (cosn(</> — the conventional 
correlation method. 

Figure || demonstrates that the results from the scalar 
product method are indeed very close to the ones of the 
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FIG. 5: The upper panel shows V2 versus centrality using the 
conventional method, where the circles and triangles represent 
t)2 with and without pt weighting, respectively. The statistical 
error is smaller than the symbol size. The lower panel shows 
the statistical error on V2 with p t weighting divided by the 
same without weighting. 



FIG. 6: The upper panel presents V2 versus centrality from 
the scalar product method (triangles) and the conventional 
random subevent method (circles). All statistical errors are 
smaller than the symbol size. The statistical error for the 
scalar product method divided by that for the conventional 
method is shown in the lower panel. 



conventional method. In this calculation, the subevents 
are generated by random partitioning. However, the de- 
tailed comparison of two results reveals a small system- 
atic difference. The difference might have origin in the 
approximations (the Central Limit Theorem) used in the 
conventional method and that are not required in the 
scalar product method. In addition, the scalar product 
method has the benefit of smaller statistical errors and 
is very simple to implement. 



III. DISTRIBUTION IN THE MAGNITUDE OF 
THE FLOW VECTOR 

In this section we study elliptic flow by analysis of 
the distribution in the magnitude of the flow vector. 
The method was used by the E877 Collaboration at the 
AGS for the first observation of anisotropic flow at ultra- 
relativistic nuclear collisions [^7|. This method is based 
on the observation that anisotropic flow strongly mod- 
ifies the distribution of the magnitude of the flow vec- 
tor J?], ||, |2(j, . Very strong flow leads to the distribu- 
tion, dP I (Q n dQn) with a local minimum at Q = 0, which 
reflects the fact that for the case of strong flow all parti- 
cle momentum unit vectors are aligned in the flow direc- 
tion. On the other hand, the non-flow effects, two and 



few particle azimuthal correlations lead to an increase in 
the statistical fluctuation width of the distribution. The 
effect can be understood by considering the flow vector 
composed of many clusters but randomly distributed in 
the azimuthal space. In the limit of large multiplicity 
and neglecting the contribution from higher harmonics 
(for a more accurate consideration se e [po[ ^(| ^9|) the 
distribution can be described by 0, |, |§T 



viM 



dP 

q n dq n 



( ln 



2o* 



q n v r 



(8) 



where Iq is the modified Bessel function. We have in- 
troduced the variable q n = Q n /VM, which greatly re- 
duces the effect on the shape of the distribution from 
averaging over events with different multiplicities. In 
a more general case using weights, one should use 
qn — Qn/(yM(wf )). In this way the width of the 
q— distribution is independent of multiplicity: 



a^= 0.5(1 + g n ), 



(9) 



with g n reflecting the change in the width of the distri- 
bution due to non-flow effects (and to some extent to the 
averaging over events with different multiplicities) . 
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FIG. 7: Reduced flow vector distributions for centrality bin 
5 plotted in two ways. Solid lines correspond to the fit with 
two parameters, Vi and gi, and dashed lines correspond to 
the fit assuming zero real flow. 



We have fitted distributions of qi , the second harmonic 
reduced flow vector, in two different ways. First, the dis- 
tributions in all different centrality bins have been fitted 
with two independent parameters, i>2 and gi- The non- 
flow contribution parameter, g2, has been found to be 
in the range of 0.18 -0.32 for all centralities except the 
most peripheral one. One should not expect a good fit 
for the most peripheral bin, for it is a mixture of events 
in a wide multiplicity range from 20 to 100. Better fit 
results for this bin could be achieved if the bin would be 
split into several sub-bins with smaller relative multiplic- 
ity variations. The relative multiplicity variation in the 
other bins is much smaller. The q distribution for the 
centrality bin 5 is presented in Fig. [?]. The two fit func- 
tions correspond to the case of a fit with two parameters, 
V2 and c/2, and to the case of a one parameter fit of g2 
for w 2 = 0. Note that the dashed curves are systemati- 
cally higher or lower than the data points in different q 
regions. In the lower part of Fig. [j] one can see that the 
anisotropic flow pushes the q distribution out to larger 
values. If the flow were great enough one could select 
events based on the q values. 

In the second method we fit q distributions in central- 
ity bins 2 to 8 simultaneously with different V2 values 
for each centrality bin but the same value of g. (This as- 
sumption is similar to the assumption of S = const in the 
previous section. See also the discussion in p8|). 
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FIG. 8: Elliptic flow as determined from the fits to the q dis- 
tributions in different centrality bins. The circles are from the 
standard method with random subevents. For the squares, all 
the centralities were fit separately. For the triangles, central- 
ity bins 2 to 8 were fit with the same value of the non-flow 
parameter. 



We find g = 0.29 ± 0.02. The results of the fits are pre- 
sented in Fig. H The deviation from the standard method 
results are due to the non-flow contributions. 



IV. FOUR-PARTICLE CORRELATIONS 

A. Motivation for Cumulants 

In experiments, it is necessary to rely on correlations 
between particles to determine the event plane since the 
reaction plane is not a direct observable. The assumption 
underlying conventional pair correlation analyses (includ- 
ing the scalar product method discussed in section II. C 
above) is that non-flow correlations of the type men- 
tioned in section I are negligible compared to the flow, 
or at most, are comparable to other systematic uncer- 
tainties. In past studies || |3(], pl| , non-flow correlations 
have been discussed with specific reference to their origin, 
such as momentum conservation, Bose-Einstein correla- 
tions, Coulomb effect, jets, resonance decays, etc. In the 
first two studies of elliptic flow in STAR (|, [lTj , the non- 
flow effect from jets and resonances was estimated using 
the approach explained in section II. A above, and this 
established an upper limit on the non-flow contribution 
to the reported Vi signal. This limit played a role in 
determining the systematic error on the published mea- 
surements. 

Anisotropic flow is a genuine multiparticlc phe- 
nomenon, which justifies use of the term collective Bow. 
It means that if one considers many-particle correlations 
instead of just two-particle correlations, the relative con- 
tribution of non-flow effects (due to few particle clusters) 
should decrease. Considering many-particle correlations, 
one has to subtract the contribution from correlations 



8 



due to lower-order multiplcts. Formally, one should use 



cumulants |2lj, |32|, |33j, 34 instead of simple correlation 
functions. Let us explain this with an example for four- 
particle correlations. The correlation between two parti- 
cles is 



{ u n,lU n2 ) 



(10) 



where n is the harmonic, and the average is taken over 
all pairs of particles in a given rapidity and transverse 
momentum region, and over all events in an event sample. 
The S n represents the contribution to the pair correlation 
from non-flow effects. Correlating four particles, one gets 



{u n ,iu n , 2 u n3 u ri 



2 • 2 



281 



(11) 



In this expression, two factors of "2" in front of the mid- 
dle term correspond to the two ways of pairing (1,3)(2,4) 
and (1,4) (2,3) and account for the possibility to have non- 
flow effects in the first pair and flow correlations in the 
second pair and vice versa. The factor "2" in front of 
the last term is due to the two ways of pairing. The pure 
four-particle non-flow correlation is omitted from this ex- 
pression — see the discussion below about the possible 
magnitude of such a contribution. What is remarkable is 
that if one subtracts from the expression ([ll]) twice the 
square of the expression ([!(]), one is left with only the 
flow contributions 

((Un,lU n ^U* n3 U* nA )) 

= (Un,l U ™,2<, 3 <,4> - 2 ( u «,l<,2> 2 = ~ v t , ( 12 ) 

where the notation ((...)) is used for the cumulant. The 
cumulant of order two is just ((it^iw* 2 )) = (u n> iu^ 2 ). 

In flow analysis, one is interested not only in so-called 
"global" flow values, but also in differential flow as func- 
tion of rapidity and transverse momentum. In a four- 
particle correlation approach, this also can be done in a 
similar manner, now correlating a particle, for example in 
a particular p t bin, with three particles from a common 
"pool". Assuming that the particle "b" is the one from 
a particular bin, one gets for a differential flow study 



Vn-bVn 



(13) 



where we have introduced the notation v n -,b for the flow 
value corresponding to the bin under study, and 6 n] b for 
the corresponding non-flow contribution. Then for the 
correlation with three particles from the pool, 

(Un,bU n ,lU* n:2 U* n3 ) = 

Vn;bvl + 2 • V^Sn-b + 2 • V n V n]b S„ + 2S n S n]b . (14) 

In this case, in order to remove the non-flow contribu- 
tion, one has to subtract from ( H ) twice the product of 
expressions ( |l0| ) and (|l3|). 

(U„,5tin,l< j2 <,3> - 2 ( u n,6<,l)<Wn,l< >2 ) = (15) 
-vh) n ;b ■ 



Assuming that the average flow value for the particles in 
the pool is known, one gets the desired differential flow 
value for the particular bin under study. 

In Eq. (pi]), we have neglected the contribution from 
the pure four-particle correlations due to non-flow effects. 
Let us now estimate the upper limit for such a contribu- 
tion. Assume that all particles are produced via four- 
particle clusters. All daughters of the decay of such a 
cluster could in principle be within 1-2 units of rapidity 
from each other. Then the contribution would be 



6//M 3 



(16) 



where M is the total multiplicity within those 1-2 units 
of rapidity, and / is ((cos2(</>i — 02)) 2 ) averaged over all 
cluster decay products. Assuming a perfect alignment, 
/ = 1, and multiplicity M — 1000, this would give us a 
possible error in v 2 measurements of the order of 



Sv 



(^ + 6/1000 3 ) 1 / 4 



V2 ■ 



(17) 



This would give only 3% relative error on v 2 signal of 
0.015, and would drops very rapidly with increasing real 
V2 signal. This calculation is for the case of 100% of 
the particle production via four-particle clusters and a 
perfect alignment of decay products. A more realistic 
scenario would give a much smaller estimate. 



B. Four-subevent method 

In order to apply the four-particle correlation approach 
to the analysis of real data, one should perform an aver- 
age over all possible quadruplets of particles in a given 
event. Bearing in mind that the average multiplicity in 
a central STAR event is well beyond a thousand, it be- 
comes a nontrivial task. The simplest solution to the 
problem is the four-subevent method where one parti- 
tions all tracks (for example, randomly) into four groups 
(subevents) and calculates a flow vector for each of the 
groups, 



Qn = ^2 '' 



(18) 



where the sum is over all particles in the group. Us- 
ing these subevents, the problem becomes much simpler 
computationally. For example, 

(«n,l«n,2< 3< 4) = (Qn.lQnjQn 3 Qn J{M 1 M 2 M Z M A )) , 

(19) 

where Mi are the corresponding subevent multiplicities. 
The cumulant calculation is straightforward: 

«Un,lUn,2<, 3 <,4» = 
,Qn,lQn,2Qn,3Qn,4: v 



M1M2M3M4 ' 

-2((g„, 1 Q; !2 /(A/ 1 M 2 )» : 



(20) 



!) 



The four-subevent method is very simple, both in logic 
and in implementation. The price for these benefits is 
lower statistical power, because the method does not take 
into account all possible quadruplets. Some improvement 
could be reached by splitting the event into more than 
four subevents and correlating all possible combinations 
of four. In the analysis of the STAR data we use eight 
subevents. A more general cumulant formalism, based 
on the cumulant generating function |3lJ offers ad- 
vantages for a four-particle analysis in the context of the 
present limited sample size. 



C. Cumulant generating function 

The cumulant and generating function approach offers 
a formal and convenient way to study flow and non-flow 
contributions systematically. Following the method of 
Rcf. the cumulant to order four is defined by 



((wiU 2 W3 W 4>> 

(uxul)(u 2 ul) 



{uiu 2 u 3 u 4 ) - 



(21) 



where, as above, the double angle bracket notation rep- 
resents the cumulant expression shown explicitly on the 
right-hand side. The subscript for the harmonic order, n, 
has been dropped. The cumulant ((uiU2U^ul)) involves 
only pure four-particle correlations, since the two-particle 
only correlations among the quadruplets have been ex- 
plicitly subtracted away. 

In the presence of flow, the cumulant becomes 

««i«a«S«4» = -«n + 0(^3+^), (22) 

where M is the multiplicity of the events, the term of or- 
der 1/M 3 represents the remaining four-particle non-flow 
effects, and the term of order v\ n jM 2 is the contribution 
of the 2n higher harmonic. The cumulant to higher or- 
ders and the corresponding generalization has also been 
determined Likewise, the cumulant of order two re- 
duces to the equivalent of a pair correlation analysis of 
the conventional type. Statistical uncertainties associ- 
ated with a cumulant analysis increase with increasing 
order from 2 to 4. 

The definition of the cumulant is simple, but it is 
tedious to calculate the moments term-by-term on the 
right-hand side of Eq. (pl|). Fortunately, the cumulant 
can be computed more easily from the generating func- 
tion 61, 



G n (z) 



M 

n 

3=1 



M 



(23) 



where z = \z\e la is an arbitrary complex number, z* 
denotes its complex conjugate. The generating function 
itself has no direct physical meaning, but the coefficients 



of the expansion of (G n ) in powers of z, z* yield the cor- 
relations of interest: 



l + (^-^}M>l"2} 

(M-l)(M-2)(M-3) , 



(24) 



)\z\ {uiu^ulul) + 



x 4M 3 

One can use these correlations to construct the cumu- 
lants. In the limit of large M, (G n ) can be used to obtain 
the cumulant generating function directly: 



M ■ 

2 k 



((Gn(z)) 1 



/At 



-1 = 



V ( fc! ) 2 



{{u l ..M k U* k+1 ...U* 2k )) . 



(25) 



The left-hand side of Eq. ( p5| ) is what is measured, and 
in order to extract the cumulants on the right, k equa- 
tions of the form of Eq. (|25|) are needed to solve for k 
undetermined parameters. This can be accomplished by 
repeating the process with k different values of \z\. It 
is found that suggested magnitudes of \z\ in Ref. pl| , 
namely roy/p with tq — 1.5 and p — 1, • • • k, are fairly 
good, since results from optimized values |35| l of r$ show 
almost no difference. Results in this paper are by default 
calculated with tq = 1.5. Since M fluctuates from one 
event to the other, for events within a multiplicity bin, 
we use the average value (M) in Eq. |25| instead of M. 

For experimental analysis, it is sufficient to take 
the first three terms in Eq. (|25|). Once the cumu- 
lant has been computed, extracting the integrated flow 
value is straightforward because, for instance, v\ — 
-((uiu 2 ulul)). 

When a non-unit weight is used, the integrated flow 
value described above becomes (w cos n<p') , which is not 
exactly v n but an approximation. However, the differen- 
tial flow can be calculated exactly (see below) no matter 
what weight is used. The integrated flow with non-unit 
weight can be obtained by integrating the differential 
flow. All integrated flow results in this paper (except 
for results from the four-subevent method) are obtained 
by integrating over the differential flow. 

For differential flow (flow in a bin of rj and/or p t ), 
Eq. ( p5| ) is replaced by 

(u d G n (z)) _ 

(Gn(z)) 



E 

k.l 



kill 



i{u d u 1 ...u k u* k+1 ...u* k+l )) 



(26) 



where Ud is the unit vector for a particle in the selected 
bin. Following a similar procedure as in the case of the 
integrated flow, the cumulant ((udUiu^u^)) is computed, 
but it now contains the angle of the one particle of in- 
terest and three other particles from the pool. Then the 
differential flow is |2l| 



{{udUiu* 2 up) 

H( U1U2M K))) 3/4 ' 



(27) 
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FIG. 9: Reconstructed V2 from the conventional method (cir- 
cles), from the 2nd-order cumulant method (triangles), and 
from the 4th-order cumulant method (stars), for simulated 
events as a function of number of embedded back-to-back 
track pairs. The horizontal dashed line marks the level of 
the true elliptic flow V2 = 0.10, as imposed on the simulated 
events, including the back-to-back track pairs. The statistical 
error is smaller than the symbol size. The multiplicity for all 
events is 500. 
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FIG. 10: Reconstructed vi from the conventional method 
(circles), from the 2nd-order cumulant method (triangles), 
and from the 4th-order cumulant method (stars), for simu- 
lated events as a function of number of embedded back-to- 
back track pairs. Unlike in the previous figure, the embedded 
back-to-back pairs are randomly distributed relative to the 
event plane, and so the true resultant V2, indicated by the 
dashed line, decreases as more pairs are embedded. The mul- 
tiplicity for all events is 500. 



Eq. (27) is for unit weight. It can be easily generalized 



for non-unit weight, and the formula still holds. 

Some detectors have substantial asymmetry in their 
response as a function of azimuth in detector coordi- 
nates, in which case it is necessary to prevent distortion 
of the measured flow signals by employing one of two 
possible compensation methods || — applying a shift- 
ing transformation which recenters Q: (<3n sinn , J , „) = 
and (Q n cosn^ n ) = (see Eq. ( g)), or applying weight- 
ing factors to force a flat \& distribution. In the present 
study, no noticeable difference is observed with and with- 
out explicit compensation for detector asymmetry, as ex- 
pected in light of the excellent azimuthal symmetry of 
the STAR TPC. All plots in this paper are made with- 
out compensation for detector asymmetry. However, it 
should be noted that cumulants, as defined by the gener- 
ating function, also correct for small anisotropies in the 
detector acceptance. For instance, the cumulant 



(28) 



amounts to an implementation of the shifting compensa- 
tion method mentioned above. 



D. Simulations 

In order to test the cumulant method as well as the 
analysis procedure, the MEVSIM (3(| event generator 
has been used to make events with various mixtures of 
flow and non-flow effects. In all cases, the number of sim- 
ulated events in a data set is 20k, and the multiplicity 
is 500. Fig. H shows one such set of simulations. Nine 
data sets with V2 = 0.10 were produced, then a simple 



non-flow effect consisting of embedded back-to-back track 
pairs was introduced at various levels, ranging from zero 
up to 80 pairs per simulated event. These pairs simulate 
resonances which decay to two daughters with a large en- 
ergy release. In Fig. ^, we consider the scenario where the 
embedded pairs themselves are correlated with the event 
plane with the same v-i = 0.10. Fig. ^ shows that the 
4th-order cumulant V2 always reconstructs the expected 
10% V2, while the V2 from the pair correlation analysis 
methods can only recover the correct input if non-flow 
pairs are not embedded. 

If back-to-back pairs are instead randomly distributed 
in azimuth, the true flow should decreases and the ex- 
pected variation can be computed knowing the number 
of random tracks. Fig. 10 shows such a simulation, and 



again it is found that only the 4th-order cumulant V2 
agrees with the expected elliptic flow, while the inferred 
V2 based on pair correlation analyses is distorted in the 
presence of the simulated non-flow effects. The role of 
resonances produced in real collisions may be closer to 
one or the other of the above two simulated scenarios, 
but in either case, the non-flow effect is removed by the 
4th-order cumulant analysis. 

In Fig. [yj consideration is given to the possible effect 
of resonances which decay with smaller energy release, 
having an azimuthal opening angle $ in the laboratory. 
The simulated events were generated with an imposed 
flow i>2 = 0.08, while in each event 50 pairs with the 
same i> were embedded, each such pair having a random 
orientation relative to the event plane. Ten data sets 
were produced, with <E> (the abscissa in Fig. [ll]) varying 
in 20° steps between zero and 180°. Again, only the 4th- 
order cumulant V2 (stars) recovers the true elliptic flow 
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FIG. 11: Elliptic flow from the conventional method (circles), 
from the 2nd-order cumulant method (triangles), and from 
the 4th-order cumulant method (stars), for simulated events 
as a function of azimuthal angle between the two tracks in 
each of 50 embedded pairs per event, with the 50 pairs each 
having random orientation relative to the event plane. The 
horizontal dashed line marks the level of the true elliptic flow 
v 2 = 0.08. 
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FIG. 12: Elliptic flow from the conventional method (cir- 
cles), and from the 2nd-order (triangles), 4th-order (stars), 
and 6th-order (crosses) cumulant methods. This is for sim- 
ulated events as a function of azimuthal angle between two 
back-to-back track pairs. The dashed line marks the level of 
the true elliptic flow. 

In order to test how the various methods respond to 
non-flow correlations associated with four-particle clus- 
ters, the simulated events in Fig. [l2| were generated with 
an imposed flow v-i = 0.10, after which 25 four-particle 
clusters were embedded in each event. Each cluster con- 
sists of two back-to-back pairs with an azimuthal opening 
angle <& between them. Seven data sets were produced, 
with <E> (the abscissa in Fig. |l^) varying in 15° steps be- 
tween zero and 90° . The clusters were oriented such that 
a track bisecting $ would contribute to the overall flow 



with V2 = 0.10. The 4th-order cumulant (stars) and the 
6th-order cumulant (crosses) both reconstruct the true 
elliptic flow (dotted line). Note that the four-particle 
correlation introduced by the clusters is 1 /M 2 times the 
pair correlation part, resulting in little difference between 
V2 from the 4th- and 6th-order cumulant methods. This 
result further illustrates the point (see also the end of 
section III. A) that non-flow effects are believed to con- 
tribute at a negligible level to the four-particle correla- 
tion, and for this reason, there may be little advantage 
in extending cumulant analyses to orders higher than 4. 
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FIG. 13: Measured elliptic flow versus centrality for Au + 
Au at ^/sjvjv = 130 GeV. The circles show the conventional V2 
with estimated systematic uncertainty due to non-flow | p7[ , 
the stars show the 4th-order cumulant V2 from the generating 
function, the crosses show the conventional vi from quarter- 
events, and the squares show the 4th-order cumulant vi from 
the four-subevent method. 



E. Results from STAR 

Figure [l3] shows measured elliptic flow versus central- 
ity, where the latter is characterized by charged parti- 
cle multiplicity n c h divided by the maximum observed 
charged particle multiplicity, n max . The conventional v% 
(circles), the 4th-order cumulant V2 from the generating 
function (stars), and the 4th-order cumulant V2 from the 
four-subevent method (squares) are compared. The cross 
symbols in Fig. [l^ represent the conventional V2 signal 
for the case where each observed event is partitioned into 
four quarter-events, which are then analyzed like inde- 
pendent events. All tracks in each quarter-event have the 
same sign of charge, and the same sign of pseudorapidity. 
Furthermore, the event plane for quarter-events is con- 
structed using only tracks with p t < 0.5 GeV/c, which 
serves to minimize the influence of non-flow associated 
with high-p t particles. It is clear that the non-flow effect 
is present at all centralities, and its relative magnitude is 
least at intermediate multiplicities. 

Figure |l4| shows i>2 as a function of pseudorapidity and 
Fig. 05] shows v 2 as a function of transverse momentum. 
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FIG. 14: Reconstructed V2 versus pseudorapidity from the conventional method (circles), from the 2nd-order cumulant method 
(triangles), and from the 4th-order cumulant method (stars), in eight centrality bins. The upper left panel shows the most 
peripheral events, and the lower right the most central. 



The eight panels correspond to the eight bins of rela- 
tive multiplicity in Fig. O but the centrality is now de- 
fined in terms of the total geometric cross section (see 
first three columns of table ||). These results illustrate 
the main disadvantage of the higher-order cumulant ap- 



proach compared with any of the two-particle methods, 
namely, larger statistical errors, and this can be seen 
to be a serious shortcoming in cases where simultane- 
ous binning in several variables results in small sample 
sizes. However, Fig. [l3] demonstrates that, especially for 
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the more peripheral bins, the statistical uncertainties for 
the fourth-order cumulant method are smaller than the 
systematic uncertainties for the two-particle methods. 

Figures |l6| and 17 are again plots of elliptic flow ver- 
sus pseudorapidity and versus transverse momentum, re- 
spectively. Here the V2 is integrated over centrality bins 2 
through 7. Bins 1 and 8 are not included in this average, 
otherwise they would significantly increase the statistical 
error on the result. The 4th-order cumulant v-i is sys- 
tematically about 15% lower than the conventional pair 
and cumulant pair calculations, indicating that non-flow 
effects contribute to v-i analyses of the latter kind. The 
vi signal based on quarter-events (as defined in the dis- 
cussion of Fig. [[3]) is closer to the 4th-order cumulant, 
although still larger on average, implying that this pair 
analysis prescription is effective in removing some, but 
not all, non-flow effects. 

Figur e O verifies that the V2(pt) curve flattens above 2 
GeV/c []3 7f1 - There is theoretical interest in the question 
of whether or not V2 (pt ) continues flat at higher p t or 
eventually goes down |38] — this issue is the subject of a 
separate analysis j3i| , and the statistics of year-one data 
from STAR is not suited for addressing this question via 
a four-particle cumulant analysis. 

Figure [l^ presents the pt-dependence of the correction 
factor for non-flow. Within errors, the relative non-flow 
effect is seen to be about the same or increasing very 
weakly from low pt through p t ~ 4 GeV/c — a some- 
what surprising result, given the presumption that the 
processes responsible for non-flow are different at low 
and high p t . Fig. 19, which presents vi from quarter- 
events divided by the conventional V2, both based on 
event planes constructed from particles with p t < 0.5 
GeV/c, offers a useful insight regarding the approximate 
Pt-independence of non-flow. This ratio roughly charac- 
terizes the contribution to non-flow from resonance de- 
cays and from other sources which primarily affect t>2 at 
lower p t , whereas non-flow from (mini)jets ought to be 
about equally present in the numerator and the denomi- 
nator of the ordinate in Fig. |l^. A comparison of Figs. [l8| 
and ^ accordingly does not contradict the implicit as- 
sumption that different phenomena dominate non-flow in 
different pt regions, and implies that the total resultant 
non-flow correction by coincidence happens to be roughly 
the same throughout the p t range under study. 

Following the approach of Section II. B, the options of 
weighting each track by either unity or p t have been com- 
pared in the 4th-order cumulant analysis. Fig. ^o] demon- 
strates that the STAR results are consistent in the two 
cases, and the pt weighting yields smaller statistical er- 
rors. All STAR results presented in this paper are com- 
puted with p t weighting unless otherwise stated. 



V. ELLIPTIC FLOW FLUCTUATIONS 

High precision results presented in this publication be- 
come sensitive to another effect usually neglected in flow 



analysis, namely, event-by-event flow fluctuations. The 
latter can have two different origins: "real" flow fluc- 
tuations — fluctuations at fixed impact parameter and 
fixed multiplicity (see, for example |Io||) — and impact 
parameter variations among events from the same cen- 
trality bin in a case where flow does not fluctuate at 
fixed impact parameter. These effects, in principle, are 
present in any kind of analysis, including the "standard" 
one based on pair correlations. The reason is that any 
flow measurements are based on correlations between 
particles, and these very correlations are sensitive only 
to certain moments of the distribution in 1*2 . In the pair 
correlation approach with the reaction plane determined 
from the second harmonic, the correlations are propor- 
tional to v 2 . After averaging over many events, one ob- 
tains (v 2 ), which in general is not equal to (v) 2 . The 
4-particle cumulant method involves the difference be- 
tween 4-particle correlations and (twice) the square of 
the 2-particle correlations. In this paper, we assume that 
this difference comes from correlations in the non-flow 
category. Note, however, that in principle this difference 



,,2\2 



7^ 0) could be due to flow fluctuations. Let 



us consider an example where the distribution in v is flat 
from v = to v = w max . Then, a simple calculation would 
lead to the ratio of the flow values from the standard 2- 
particle correlation method and 4-particle cumulants as 
large as {v 2 ) 1 ' 2 /(2{v 2 ) 2 - {v 4 )) 1 ' 4 = h 1 ' 4 r> 1.5. 

In this study, we consider the possible bias in elliptic 
flow measurements under the influence of impact param- 
eter fluctuations within the studied centrality bins. The 
largest effect is expected within the bin of highest mul- 
tiplicity, where the impact parameter and V2 are both 
known a priori to fluctuate down to zero in the limit of 
the most central collisions. These fluctuations lead to 
bin-width-dependent bias in the extracted V2 measure- 
ments. 

In section III, two approximations were made in order 
to extract the final flow result, 
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and 



Taking into account the centrality binning fluctuation on 
flow, namely a 2 2 and <j 2 n , 



\V„. = <J„2 



and 



,,2 \ 



(Vn) 2 



and Eq. <X2W) becomes 



2a 2 vl 



2 



(29) 



which is a function of v n and is solvable for v n , if a% and 
a 2 ? are known. A method of calculating both a 2 n and 
a 2 ? is now presented. 

First, we need to parameterize function of 

impact parameter, b. Consider a polynomial fit v n — 
ciq + a\b + ... + a^b 6 , in which case the measured flow is 
(v n ) — a + ai(b) + ... + a 6 (6 6 ). The various averages (6), 
(6 2 ),... (6 12 ) can be estimated in each centrality bin from 
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FIG. 15: Reconstructed V2 versus pt from the conventional method (circles), from the 2nd-order cumulant method (triangles), 
and from the 4th-order cumulant method (stars), in eight centrality bins. The upper left panel shows the most peripheral 
events, and the lower right the most central. 



filtered HIJING events. The parameters a* have been 
determined by minimizing \ 2 m a fit to the eight i>2(?ich) 
measurements. In addition, the fit is constrained to go 
through v 2 — at b — and at b max = 14.7 fm Q. The 
variation of b max within ±0.5 fm has a negligible effect 



on V2(b) at b < 12 fm. Fig. pi] shows the resulting curve: 

v 2 (b) = -0.000394 b + 0.00210 b 2 
-0.0000706 b 3 - 0.0000320 b A + 
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FIG. 16: Elliptic flow versus pseudorapidity from the conven- 
tional method (circles) , from the 2nd-order cumulant method 
(triangles), from quarter-events (crosses), and from the 4th- 
order cumulant method (stars), averaged over all centralities 
from bin 2 through 7, as defined in Figs. and 



FIG. 18: The ratio of V2 from the 4th-order cumulant di- 
vided by «2 from the conventional method as a function of pt , 
averaged over all centralities from bin 2 through 7, as defined 
in Figs. and [THJ. 




3.5 4 4.5 5 

p, (GeV/c) 

FIG. 17: Elliptic flow versus transverse momentum from 
the conventional method (circles), from the 2nd-order cumu- 
lant method (triangles), from quarter-events (crosses), and 
from the 4th-order cumulant method (stars), averaged over 
all centralities from bin 2 through 7, as defined in Figs. |l4| 
and|l|. 
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FIG. 19: The ratio of V2 from quarter-events divided by the 
conventional V2 as a function of pt- In both cases, event planes 
were constructed from low pt (< 0.5 GeV/c) particles. The 
data are averaged over all centralities from bin 2 through 7, 
as defined in Figs. [l4| and [lit 



0.00000358 b 5 - 1.174 x 10~ 7 b 6 , (30) 

where it is assumed that b is in fm. In principle, the final 
corrected ^(^ch) should be determined iteratively, but 
the result is stable on the first iteration. 
Next we consider 

°u„ = ( V l) - M 2 = 

{a 2 + 2a a 1 (b) + ...a 2 e (b 12 ))- 

(ao + ai <6) + ...a 6 (6 6 )) 2 , (31) 

and again the various averages of powers of b can be 
estimated using HIJING. 

After computing a 2 n , a 2 2 , and obtaining ?; moas from 
the four-particle correlation method, Eq. (|2^) can be 



solved to extract the v n corrected for impact parameter 
fluctuations. The vi bias is found to be entirely negligi- 
ble in all the studied centrality bins except for the most 
central, where the correction is about a factor of two (see 
the leftmost bin in Fig. |2l]). In the present analysis, even 
a factor of two is not significant due to the large statis- 
tical error on V2 for maximum centrality. However, the 
correction to Vi resulting from finite centrality bin width 
at maximum centrality has been determined with lower 
uncertainty than V2 itself, and will become important in 
future studies with large samples of events. 

Real event-by-event fluctuation in the flow coefficients 
would also make the four-particle values lower than the 
two-particle values. At the moment, there is no way to 
calculate this effect, although it is expected to be small. 
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{rich} 


(j^ch f Umax } 


cross section 


(6) (fm) 


RMS (6) 


(e) 


(V2) 


5333 


0.060 


53 - 77% 


12.23 


0.99 


0.420 


0.052 ±0.012 


14038 


0.160 


41 - 53% 


10.36 


0.70 


0.415 


0.055 ±0.003 


2271? 


0.258 


31 - 41% 


9.06 


0.68 


0.371 


0.053 ±0.001 


319S 


0.363 


24 - 31% 


7.91 


0.64 


0.319 


0.051 ±0.001 


415fl 


0.472 


16 - 24% 


6.80 


0.70 


0.261 


0.043 ± 0.001 


519g 


0.590 


10 - 16% 


5.56 


0.72 


0.197 


0.035 ± 0.002 


622i 


0.708 


5 - 10% 


4.26 


0.80 


0.131 


0.023 ± 0.002 


746i 4 


0.849 


top 5% 


2.53 


1.00 


0.058 


0.012 ±0.015 



TABLE I: Tabulated values of observed charged particle multiplicity, n c h/n max , centrality in terms of percent of total geometric 
cross section, impact parameter with spread (root mean square) inferred from HIJING, the initial spatial anisotropy e, and the 
final corrected elliptic flow based on 4th-order cumulants. 
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FIG. 21: V2 as a function of impact parameter. The data 
points are shown at the values of (b) for a given centrality 
bin. For the description of the fit procedure see text. The tri- 
angles are the final 4th-order cumulant data after correction 
for fluctuations as described in Section IV, while the circles 
show the 4th-order cumulant data before this correction. The 
dashed lines represent the estimated uncertainty in the pa- 
rameterization represented by the solid curve. 



n cl/ n max 



FIG. 20: The upper panel shows measured vi from 4th- 
order cumulants versus centrality with pt weighting (stars) 
and unit weighting (circles) . The bottom panel is the ratio of 
error from pt weighted V2 to that of unit weighted V2 ■ 



VI. THE CENTRALITY DEPENDENCE OF 
ELLIPTIC FLOW 

The centrality dependence of elliptic flow is a good 
indicator of the degree of equilibration reached in the 
reaction [13). Following Ref. [fy], we compute the 
initial spatial eccentricity for a Woods-Saxon distribution 
with a wounded nucleon model from 

e = (y 2 ) (^ 2 ) 
(y 2 ) + (x 2 ) 



where x and y are coordinates in the plane perpendicular 
to the beam and x denotes the in-plane direction. The 
method of calculation of epsilon is the same as used for 
the hydro values [Q . The ratio v% /e is of interest because 
it has been argued to be independent of centrality in a 
hydrodynamic model with constant speed of sound 
In hydrodynamic model calculations using an equation 
of state with a phase transition (sound speed is not con- 
stant) this ratio does change as a function of centrality, 
however within the 10% level [Q. Hydrodynamics rep- 
resents one possible limiting case in describing nuclear 
collisions — the limit where the mean free path for inter- 
action of the constituents represented by the fluid cells is 
very small compared with the region of nuclear overlap. 
The opposite limit, where the mean free path is long (or 
at least comparable to the dimensions of the nuclear over- 
lap region) is normally known as the Low Density Limit 
(LDL). In nuclear transport models, the mean number 
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of hard binary interactions per particle is typically small, 
and the predictions of these models tend to be closer to 
the low density limit than the hydro limit. In order to 
judge the proximity of measured flow data to either of 
these limits, it is useful to plot, as in Fig. ^2], v-i/e versus 
charged particle density in the form (dN/dy)/ S, where 
dN/dy is rapidity density, and the area of the overlap re- 
gion is S — ~K\J (x 2 ) (y 2 ) as computed above. Since V%/e is 



proportional to (dN/dy)/ 'S in the LDL case J_8|, |42j, this 
form of plot offers meaningful insights without reference 
to detailed theoretical models. 
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FIG. 22: V2/e as a function of charged particle density in An 
+ Au collisions. Data are from E877 at the AGS (squares), 
NA49 at the SPS (circles), and STAR at RHIC (stars). The 
AGS and SPS data have been obtained by conventional flow 
analysis. The STAR measurements are at t/snn = 130 GeV, 
and correspond to the final corrected elliptic flow based on 
4th-order cumulants, and we assume dN/dy — 1.15dN/dy. 
The horizontal shaded bands indicate the hydrodynamic lim- 
its for different beam energies flit . 

Figure |2| presents Au + Au data from AGS/E877 
@, from NA49 Q, as well as the current STAR 
measurements based on 4th-order cumulants, corrected 



for fluctuations as detailed in Section IV. Alternative 
forms of the centrality dependence readily can be gen- 
erated using the tabulated quantities presented in Ta- 
ble I. Generally, the current STAR results underline 
the need for much increased statistics, particularly for 
the most central collisions. Within the uncertainties, a 
smooth trend of increasing v-i/e with increasing central- 
ity (larger (dN/dy)/S) is observed, without the obvious 
kink that has been suggested as a phase transition sig- 
nature |l8|, [hi). Another proposed phase transition sig- 
nature which is not favored by the data is a few percent 
rise in V2/S with decreasing centrality 44 . It is note- 



worthy that the v^/e values reached in the most central 
RHIC collisions are consistent with the hydrodynamic 
limit H ^J, [ItJ , whereas v-i /e in central collisions at AGS 
and SPS is significantly lower. It is also worthy of note 
that while the roughly linear relationship between v-i/s 
and (dN/dy) / S across the presented beam energies and 
centralities is consistent with the LDL picture [pL8| , the 
measured V2(pt) Fig. [l 5| c a nnot be explained by current 
LDL implementations |48|, and is much closer to hydro- 



dynamic calculations up to 2 GeV/c 48 



VII. CONCLUSION 

In this work, we provide details of the approach for 
treating non-flow correlations within the framework of 
the standard elliptic flow analysis method based on par- 
ticle pairs. We also compare the standard method with a 
new and simpler pair analysis based on the scalar product 
of flow vectors. The latter yields a 15 - 35% reduction in 
statistical errors, with the best improvement occurring 
in the case of the most central and the most peripheral 
events. 

It is concluded that four-particle correlation analyses 
can reliably separate flow and non-flow correlation sig- 
nals, and the latter account for about 15% of the ob- 
served second-harmonic azimuthal correlation in year- 
one STAR data. The cumulant approach has demon- 
strated some advantages over the previous alternatives 
for treating non-flow effects. In particular, 4th-order cu- 
mulants allows us to present vi measurements fully cor- 
rected for non-flow effects, in contrast to the earlier anal- 
yses where the non-flow contribution was partly removed 
and partly quantified by the reported systematic uncer- 
tainties. It is observed that non-flow correlations are 
present in -^sjvw = 130 GeV Au + Au events throughout 
the studied region < 1.3 and 0.1 < p t < 4.0 GeV/c, 
and are present at all centralities. The largest contribu- 
tion from non-flow correlations is found among the most 
peripheral and the most central collisions. 

On the other hand, a 4th-order cumulant analysis is 
subject to larger statistical errors than a conventional 
pair correlation analysis of the same data set. The to- 
tal uncertainty on the 4th-order analysis, including both 
statistical and systematic effects, is smaller for year-one 
STAR data except in the most central and peripheral 
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panels of Figs. |lj and Il5[ In the case of future studies of 
larger numbers of events, a higher-order analysis should 
provide an advantage in all cases. 

Fluctuations within the studied multiplicity bins have 
the potential to bias elliptic flow results. This bias has 
been estimated and found to be entirely negligible except 
for the most central multiplicity bin, where the correction 
is about a factor of two. In the present analysis, even this 
large a bias is only marginally significant, but again, this 
correction will presumably be important in future studies 
with much improved statistics. 

We present STAR data for V2 /e — elliptic flow in var- 
ious centrality bins, divided by the initial spatial eccen- 
tricity for those centralities. Mapping centrality onto a 
scale of charged particle density enables us to study a 
broad range of this quantity, from peripheral AGS colli- 
sions, through SPS, and ending with central RHIC col- 
lisions. Within errors, the STAR data follow a smooth 
trend. No evidence for a softening of the equation of state 
or for a change in degrees of freedom has been observed. 
The three experiments at widely differing beam energies 
show good agreement in V2 / £ where they overlap in their 
coverage of particle density. The pattern of t^/e being 
roughly proportional to particle density continues over 
the density range explored at RHIC, which is consistent 
with a general category of models which approximate the 
low density limit as opposed to the hydrodynamic limit. 



Nevertheless, v-ije at STAR is consistent with having just 
reached the hydrodynamic limit for the most central col- 
lisions. 
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